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1 Motivation 



' We review in brief the threshold expansion, a method to perform the expansion of Feynman 

^ , integrals near the heavy quark-antiquark threshold, and its relation to the construction of two 

effective theories, non-relativistic QCD (NRQCD) and potential-NRQCD. We then summarize 
' recent next-to-next-to- leading order results on the decay J/^ — > , the bottom quark 

^Sj , mass from T sum rules and the top-anti-top production cross section near threshold in e^e~ 

collisions. 

(N 
> 

(N 
^ , 

, In this talk we discuss systems of the form QQ + X, where Q is a heavy quark of mass m and 

' X a collection of massless particles, in the kinematic region, where the total invariant mass q'^ 

Q>^ ■ of the system is close to Am?. The physics of such systems is characterized by the fact that in 

the reference frame with q = the heavy quark velocities are small, i.e. u <C 1. 
0^1 The threshold region is evidently sensitive to the mass of the quark and physical quantities 

Q^' that probe the threshold region can therefore be used to determine the bottom and top quark 

O . mass. In case of bottom quarks the threshold region is populated by narrow Upsilon resonances 

and the cross section e^e~ bb + X cannot be predicted locally. However, dispersion relations 
equate an average over Upsilon resonances and the continuum to the derivatives of 5-<Tuark 
^ . current correlation function at = 0, which can be calculated in perturbation theory Er@. In 

case of top quarks the rapid decay of the top prohibits the formation of long-lived resonances §. 
For mt = 175 GeV a resonance- like structure remains visible in the energy dependence of the 
cross section e'^e~ tt -\- xHi. Perturbation theory can still be applied in the vicinity of the 
'pseudo-resonance', but not arbitrarily close to threshold. The shape of the tt cross section near 
threshold, to be measured at a Next Linear Collider, is believed to provide us with the most 
accurate determination of the top quark mass, if the strong interaction corrections are indeed 
well understood. 

In addition to quark masses as parameters of the standard model, there are more intrin- 
sically QCD-related problems involving heavy quarks near threshold. Heavv quarkonia are 
non-relativistic and their production and decay can be treated as an expansionQ in v. One can 
also consider heavy quark production in hadron-hadron collisions. In this case one is mainly 
concerned with the resummation of logarithms of v that arise as a consequence of soft and 
collinear gluon emission, predominantly from the massless ioitjal state particles. This last ap- 
plication will not be discussed here (see, for instance, Refs. Incidentally, we note that at 
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Tevatron centre-of-mass energies the average velocity squared of top quarks is (f^) ~ 1/2 and 
the non-relativistic dynamics discussed subsequently is probably not important there. 

A non-relativistic system involves more than one momentum scale, related to the small pa- 
rameter V. The first part of the problem consists of organizing the calculation in a systematic 
expansion in v and has been addressed in several papers recently ll2rllZI, especially in the con- 
text of diipenaipnal regularization. These works clarify the definition of non-relativistic QCD 
(NRQCD)lla~c3 with a matching prescription based on dimensional regularization and show that 
one can proceed to a second effective theory by integrating out more scales. This development 
is summarized in the first part of the talk. 

The second part summarizes recent next-to-next-to-leading order (NNLO) results on the 
phenomenological applications mentioned above. 'NNLO' in this context refers either to fixed- 
order calculations in ag required for matching effective theories, or the resummation of 
Coulomb-enhanced corrections to all orders required for the solution of the low-energy 
problem. 



2 Complications near the threshold 

The need for resummation arises, because even when is small the effective interaction between 
a heavy quark and anti-quark becomes strong at small relative velocity. The exchange of a 
Coulomb gluon - the 00-component of the propagator for a gluon with momentum of order mv 
and energy of order mv'^ - leads to an effective coupling of order Tras{mv) /v. Resummation leads 
to weakly coupled bound states (-Ebind ~ '"T'^'^ <^ m) - the Coulomb bound states analogous 
to hydrogen and positronium. In first approximation the physics is indeed exactly as in QED. 
However, due to massless quarks and gluons, the coupling runs below the scale m in QCD. This 
leads to complications, and differences in the power counting, when mv'^ ~ Aqcd- 

Defining v = {1 — Am? / q^Y^"^ , the cross section e^e" QQ + X can be expanded in a 
double series in Og and v: 

^QQ =^Y1 ^kiasV^ X logs of V. (1) 
k=0 l=-k 

A NNLO calculation in the kinematic region where ~ u has to account for all terms with 
k + I < 2. Powers of v arise from ratios of momentum scales. We have to disentangle the 
contributions from the different scales in order to be sure that for a high-order loop graph, 
which cannot be calculated exactly, we have taken into account all terms with k + I < 2. 

Eventually we will be led to calculating diagrams with Coulomb Green functions rather than 
free Green functions. There are, however, two 2-loop calculations in conventional perturbation 
theory that enter at NNLO: (a) The 2-loop correction to the Coulomb potential, i.e. 02 in 

1 k 

It has been calculated in Ref. 0. (b) The 2-loop correction to the matching of the vector heavy 
quark current in QCD to its non-relativistic analogue, i.e. C2 in 

QliQ = CfcaJ^ iP^aiX + ■■■, (3) 

where ip and x 2-spinors. Onlv the spatial component is needed at NNLO. The coefficient 
C2 has been calculated in Refs.EBcl. Note that the series expansion of the Coulomb potential 
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is uniquely specified (at least to order a^) by choosing a scheme for coupling renormalization. 
But the non-relativistic cur rent is not conserved and one has to choose a factorization scheme 
to define the Cfc. In Refs.c3IH3 the MS scheme is used. 

3 Scales and their Separation 

We now discuss how to construct the expansion (I), first at a given loop order k, then including 
resummation. 

3. 1 Threshold expansion of Feynman integrals 

Take a loop integral that contributes to the heavy quark cross section. To be precise, we consider 
the two-point function of the current @ and take the imaginary part at the end. This avoids 
calculating the threshold expansion of phase space integrals. The threshold expansion llj is 
a method to calculate the expansion in v without calculating the full integral. Inspection of 
the denominators of the Feynman integrand shows that there are four essential regions of loop 
momentum: 



hard (h) 
soft (s) 
potential (p) 
ultrasoft (us) 



Iq ^ m, I m, 
Iq mv, I ~ mv, 

Iq ~ mv'^, I ~ mv, (4) 
Iq ~ mv'^, I ~ mv'^. 



The threshold expansion is constructed by writing the diagram as a sum of terms that follow 
by dividing each loop momentum integration into these four regions. More precisely, one has to 
account for the possibility that for example li and I2 are hard, but their sum is not. That is, one 
should sum over one-particle- irreducible subgraphs. Note that the division is done implicitly, 
through the expansions of the propagators. No explicit cut-offs are needed. The soft region 
has often been omitted in the discussion of non-relativistic dynamics, probably because it gives 
rise only to instantaneous potentials as we discuss below. However, it also gives rise to the 
(standard) evolution of the strong coupling between the scales m and mv, as can be seen from 
the fact that a light quark loop inserted into a potential gluon line can only be hard or soft. 

Then, in each term, the propagators can be simplified. The propagator of a heavy quark 
with momentum {q/2 + Iq,p + I) is 
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(5) 



ll -l^ - 2p- 1 - p'^ + qlo + y + ie' 

and we will assume that p scales as mv and y = — as mv'^. When / is hard, we expand 
the terms involving p and y and the leading term in the expansion scales as v^. When I is 
soft, the term qlo is largest and the remaining ones are expanded. The propagator becomes 
static and scales as v~^. When / is potential, the propagator takes its standard non-relativistic 
form after expansion of Iq and scales as Massive particles can never be ultrasoft. The 

gluon propagator takes its usual form, when the gluon line is soft and ultrasoft and scales as 
v"'^ and v~'^, respectively. If the gluon momentum is potential, one can expand Iq and the 
interaction becomes instantaneous. If we add the scaling rules for the loop integration measure, 
d'^l ~ 1(h), f^(s), (p), v'^(us), we can immediately estimate the size of the leading term 
from a given region. Because all terms that are small in a given region are expanded in the 
Feynman integrand, each term in the resulting sum contributes only to a single power of v 
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(multiplied, in general, by powers of logarithms of v). It is important that the power behaviour 
can be determined before calculating. It is also important that the integrals that appear in each 
term of the expansion are much easier to calculate than the original integral, because they are 
homogeneous in that is, contain only a single scale. In particular, any term in the expansion 
of massive 2-loop 3-point integrals is calculable EJIHa. 

The diagrammatic rules for the expansion of propagators (and vertices) can be considered 
as following from an effective Lagrangian. Because different loop momentum regions do not 
overlap, one can introduce separate fields for hard, potential, etc. quarks and gluons. (For the 
coupling of potential quarks to ultrasoft ffluons, the threshold expansion entails the multipole 
expansion.) This has been done in Ref. Ii3 and, in part, in Refs.li^JtJ. One can then think of 
integrating out first the hard fields, then soft fields and potential gluons. 

One can use Coulomb gauge or a covariant gauge. (In covariant gauge the ghost fields are 
treated like other massless fields. The can be hard, soft, potential and ultrasoft.) The hard 
region is very inconvenient to calculate in Coulomb gauge. But Coulomb gauge makes gauge 
cancellations manifest that occur in the coupling of soft and ultrasoft gluons to heavy quarks. 
In fact, one can use different gauges for different regions with one exception: The distinction 
between soft and potential gluons is not gauge-invariant. This should be no surprise. When one 
integrates out soft gluons alone, one has to compute graphs with external potential gluons. But 
potential gluons are off their mass shell. On the other hand, when one integrates out soft and 
potential gluons together, only on-shell graphs have to be considered. 



3.2 NRQCD 

First we integrate out the hard region. This means that we assume that all loop momenta 
are hard and the external momenta soft, potential or ultrasoft. That is, we Taylor-expand 
the integrand in the external momenta of the graph. The result is obviously polynomial in 
the external momenta and can therefore he written as a local operator. This yields the non- 
relativistic QCD (NRQCD) Lagrangian IiUhoI 

- ^ '4)Ud ■ E - E ■ D^ij - ^ '4)'^& ■ (d X E - E X D)^) 

-|- antiquark terms -|- £iight (6) 

and an expansion of QCD operators (such as the vector current) in terms of non-relativistic 
fields. We have written down only those terms of the Lagrangian that are needed for the NNLO 
calculations described below, provided one treats the difference between Q;s(m) and as{mv'^) as 
small. Otherwise one should use the renormalization group to sum up the leading logarithms 
sensitive to the scale mv^ in the coefficient functions. 

Note that the threshold expansion provides a matching prescription, in which one does not 
have to calculate NRQCD diagrams explicitly. The QCD diagram is already broken up into 
the contributions from different scales and the hard regions are exactly the relativistic eff'ects 
which are contained in the coefficient functions of the effective Lagrangian. The matching 
prescription is very simple. Despite the fact that the QCD diagrams are divergent, when the 
relative momentum of the heavy quarks goes to zero, the matching coefficients are given by 
the Taylor expansion around zero relative momentum. For the single-heavy quark sector, this 
matching prescription coincides with that of Ref. El. 
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Note also that perturbative calculations with the NRQCD Lagrangian have up to now been 
mainly done with a cutoff regularization, either in QED or lattice QCD, while we assume di- 
mensional regularization here. The change is not quite trivial. For, if we calculated a NRQCD 
integral according to the Feynman rules of its Lagrangian, we would obtain an incorrect result, 
because dimensional regularization treats the UV cutoff of NRQCD as larger than m. We can 
use dimensional regularization provided we expand the integrand before integration. The pre- 
cise prescription is again specified by the threshold expansion. The NRQCD integral is given 
by expanding the integrand according to the rules for the soft, potential and ultrasoft region. 
Indeed, since the hard region is accounted for by the coefficient functions, this reproduces the 
original QCD diagram. 

The interaction terms in NRQCD do not have a definite scaling in v, because a gluon field 
can be soft, potential or ultrasoft, and a quark field can be soft or potential and different scaling 
rules apply to each of these. 

3.3 Potential NRQCD 

If the scale mv is already non-perturbative, we stop with NRQCD. If it is not, the threshold 
expansion suggests that one can integrate out the soft region together with potential gluons to ar- 
rive at another effective Lagrangian. We call the effective theory potential NRQCD (PNRQCD), 
following Refs.O'Ej, where the corresponding tree- level Lagrangian was considered first. We can- 
not integrate out potential quarks, because the Green functions relevant to our applications have 
external potential quark lines. 

To integrate out the soft contributions we consider (NRQCD) graphs in which all momenta 
are soft and the external momenta potential or ultrasoft. The graph is Taylor-expanded in the 
external ultrasoft momenta and in the zero components of external potential momenta. But 
it cannot be expanded in the spatial components of potential external momenta, because they 
are not small compared to the spatial components of the loop momenta. Hence the result is 
non-polynomial in the spatial components of the potential momenta. This is an instantaneous, 
but spatially non-local interaction. We call such interactions potentials. Potential gluons have 
propagators with Iq expanded. They also give rise to potentials. 

PNRQCD contains only potential quark fields and ultrasoft gluon (massless quark, ghost) 
fields. NRQCD is matched on PNRQCD order by order in as- Consider quark-antiquark 
scattering at small relative momentum. At leading order in Os the quark and antiquark interact 
by the exchange of a potential gluon. The leading term in v yields the Coulomb potential at 
order a^. The corresponding non-local operator is 



The corrections of order v"^ are known as the Breit potential. At order one has to compute the 
soft and potential contributions to the 1-loop NRQCD graphs and subtract the PNRQCD graphs 
constructed from the order-o^ potentials in the PNRQCD Lagrangian. The soft contributions to 
NRQCD graphs have no analogue in PNRQCD and renormalize the PNRQCD interactions. The 
1-loop correction to the Coulomb potential is generated in this way together with a potential of 
form a^/(mr^) and higher order terms in v. For the potential contributions to NRQCD graphs 
it is necessary to perform an explicit matching to avoid double counting. As mentioned above, 
the contributions from soft and potential gluons may look different in different gauges. But their 
sum and hence the PNRQCD Lagrangian is gauge- invariant. 
In general the PNRQCD Lagrangian can be written as 




(7) 



-CPNRQCD = ^'^ 



■NRQCD 



■non— local ? 



(8) 
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where >Cnon-iocai collects all non-local interactions. The local interactions are exactly those of 
NRQCD, but the interpretation is different, because only ultrasoft gluons are left over. In loop 
graphs constructed from >C^j^qqj-, the gluon propagators are always expanded according to their 
ultrasoft scaling rule, while in loop graphs constructed from >Cnrqcd gluons are also soft and 
potential. The prime reminds us of this difference. 

Because only potential quarks and ultrasoft gluons are left in PNRQCIl, the interaction 
terms have definite scaling rulss. They agree with those given in Refs.ll°Hla. Note that the 
NRQCD scaling rules of Ref. Ej are really those of PNRQCD. A potential quark propagator 
in coordinate space scales as v^, so a quark field in PNRQCD scales as w^/^. Comparing the 
scaling of ip^d^ip with the scaling of (0), we find that the former scales as and the latter 
scales as a<jU^. As is well-known, the Coulomb interaction cannot be treated as a perturbation 
when V ~ as{mv). Note, however, that the matching on PNRQCD can be done by treating the 
Coulomb interaction as a perturbation. The unperturbed PNRQCD Lagrangian is 



^tT^V] in (-^) [x^T^x] (0) + Clll^,. (9) 



One can rewrite this in terms of a 'tensor field' [t/j that depends on the cms and 

relative coordinates. The unperturbed Lagrangian describes free propagation (with mass 2m) 
in the cms coordinate. The propagation of [ip (8> x^{Ri^) iii its relative coordinate is given by 
the Coulomb Green function of a particle with reduced mass m/2. In calculating diagrams with 
Coulomb Green functions, one sums corrections of order (as/v)" to all orders. The remaining 
terms can be treated as perturbations in v and as around the unperturbed Lagrangian. These 
are calculations familiar from QED bound state problems. What is new, in our opinion, is that 
we understand how to perform such calculations systematically in dimensional regularization, 
without double counting, and, if necessary, including retardation effects (graphs with ultrasoft 
lines) . 

When an ultrasoft gluon line with momentum I connects to a quark line with loop momentum 
k — l/2 for the incoming and k + l/2 for the outgoing quark line, the threshold expansion instructs 
us to expand the quark-gluon vertex and quark propagator in l/k ~ v. All gluon interaction 
terms in /^nrqcd should be understood as multipole-expanded, for instance 

ip^^A-dip] (x) = ^'^{t,x)A{t,0) ■ dtp{t,x) +11^^(1, x) {x-d)A{t,0) -d^it^x) + ..., (10) 

and likewise for all other interactions. 

Up to this point we have neglected the fact that in QCD - contrary to QED - the coupling 
constant evolves below the scale m. When mv"^ <C m, but as{mv^) <C 1, one can sum up 
logarithms of but otherwise the power counting remains unaffected. In particular, only the 
Coulomb interaction has to be treated non-perturbatively. When mv"^ ~ Aqcd the situation 
changes, because ultrasoft gluons couple with unit strength, since asimv^) ~ 1. The coupling 
to heavy quarks is still small, of order v at least, but the self-coupling of gluons is unsuppressed. 
An ultrasoft gluon propagator in coordinate space scales as f'*, hence the gluon kinetic term 
scales as v^. The three-gluon and four-gluon vertices also scale as v^. When as{'mv'^) ~ 1, 
ultrasoft gluons enter the calculation of Rqq as a non-perturbative contribution of relative 
order and 'retardation effects' cannot be neglected at NNLO. A perturbative treatment of 
the problem cannot be extended beyond NLO, because the unperturbed PNRQCD Lagrangian 
must be modified to contain £iight rather than /^[[g^f The unperturbed problem is then no longer 
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exactly solvable. This is why the energy spectra of charmonia and bottomonia are not exactly 
Coulomb-like. 



4 J/?/; ^ 

The two-loop short-distance correction to the leptonic decay of an S-wave quarkonium state 
such as J/ip, tp' and T{nS) has been analyzed in Ref. c3. The decay rate can be expressed as 



, fc=0 



4nelal^ 12 |M'(0)|2(/.) 



3M, 



neglecting relativistic corrections which can be added systema ticallyl at the expense of further 
non-perturbative parameters in addition to the wave-function at the origin. The coefficients Ck 
are those of (^) (up to a normalization) and ^'(0) is related to the vacuum-to- J/V' matrix element 
of the non-relativistic current. We have cq = 1, ci 
MS factorization scheme reads E3o 



(11) 




and the 2-loop coefficient in the 
C(3) 



+ CfCa<tt' t In ^ +--T 



+ CfTf 




(12) 



The calculation amounts to picking up the hard contributions only of the 2-loop three-point 
integrals contributing to c -|- c — > 7* . Numerically 



k=0 



asirric 



vr 



1 



37r 



44.55 - 0.41n/ - 25.59 In 



TJln 



+ .... 



(13) 



Take /i = nic for the factorization scale. Before squaring the coefficient function, the 1-loop 
correction is —25%, but the 2-loop correction amounts to —50%. Even for bottomonium, the 
2-loop correction is as large as the 1-loop correction. 

At two loops the wavefunction at the origin becomes factorization scale and scheme-depen- 
dent. The anomalous dimension is very large. The leptonic width is an important observable 
with respect to tuning the parameters of potential models. The large scheme-dependence is a 
problem for potential models, because it is not clear which scheme the wavefunction at the origin 
in potential models corresponds to. 

The above result suggests that the perturbative expansion is not reliable, so that pertur- 
bative factorization would not work quantitatively. But the large coefficients could also be the 
consequence of a 'bad' factorization scheme. We will know only once a second quarkonium decay 
such as rjc — > 77 is computed to second order. With present analytic methods this seems to be 
a challenging piece of work. 



5 tt Production 



5.1 The total cross section 



The total cross section for e~^e tt + X has been calculated at NNLO in Refs. rW 1 for the 
vector coupling of the tt pair. (The axial-vector contribution is a NNLO effect but has not been 



7 




Figure 1: The cross section as function of y/s — 2mt (in GeV) in LO (dotted), NLO (dashed) and NNLO 
(sohd) for three choices of renormahzation scales each. Parameters: = 175 GeV, Fj = 1.43 GeV and 

as{mz) = 0.118. Figure reprinted from Rcf. EZl. 



considered in Refs.rnrl.) This can be done by combining the hard matching coefficient ( [T^ ) with 
the integrals over Coulomb Green functions in dimensional regularization. In practice, Ref. 
used a different factorization scheme by comparing the resummed cross section with the cross 
section at the threshold at fixed order a^B. The final result is independent of this choice. 

The top quark is unstable on a scale comparable to the scale mja^ and finite width effects 
are essential in the threshold region. Following Ref. H, the finite top quark width is taken into 
account by the substitution 

E = ^/^-2mt^E + irt (14) 

in the argument of the Coulomb Green function, where E is the energy measured from the 
threshold, defined by twice the top quark pole mass. Note that E depends on the renormalization 
convention for the top quark mass and is not a physical quantity, contrary to the cms energy ^/s. 
The prescription ( p!^ ) has been justified in Ref. □ at LO. It is probably not justified at NNLO. 
In general, one can expect that the finite width inhibits the radiation of ultrasoft gluons from 
top quarks but affects the potential interactions less strongly. 

The result for the cross section in units of the point cross section as a function of E is 
reproduced in Figure |l|. We note that the NNLO calculation has an as large effect on the 
height of the peak as the NLO calculation. Furthermore both shift the location of the peak by 
somewhat less than a GeV. If the QCD corrections do not converge, this implies an uncertainty 
in mt of almost 500 MeV, with frustrating consequences for precision studies of the ti threshold 
at a Next Linear Collider. 



5.2 Which mass? 

There is reason to believe that the situation is not as bad. When one discusses uncertainties 
in quark masses in the range of a few hundred MeV, it is important to ask how sensitive a 
particular mass renormalization prescription is to long-distance QCD effects. One has to be 
particularly careful about this, if the physical process is intrinsically less long-distance sensitive 
than the mass renormalization convention one is about to use. 
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For the pole mass of a stable quark long-distance sensitivity causes an uncertainty of order 
Aqcd^S- The analysis of Refs.HEl goes through unmodified^ for an unstable quark, when the 
pole mass is defined as the real part of the complex pole position in the top quark propagator. In 
perturbative calculations this long-distance sensitivity shows up as large radiative corrections, 
when the top quark pole mass as an input parameter is fixed as done above. It is advantageous to 
fix instead a mass renormalization scheme which is less sensitive to long-distances, provided one 
can show that the large corrections cancel then. (Effectively, this amounts to comparing a LO, 
NLO etc. calculation in terms of the pole mass at somewhat different input values of the pole 
mass. Because of the definition (14) of E, this leads to a shift in the horizontal scale of Figure ||, 



which compensates the shift in the LO, NLO and NNLO curves. It seems preferred to plot 
the cross section as a function of the physical parameter y/s rather than the scheme-dependent 
parameter E.) 



Indeed the analysis of 
systematically large terms 



38h 



40 



order radiative corrections to the Coulomb potential reveals 
In part, these large corrections are a consequence of long- 
distance sensitivity. The long-distance sensitive contributions to the Coulomb potential can 
be shown to cancel exactly against those in the pole mass renormalization. Contrary to 
intuition the threshold cross section is less sensitive to long distances than the pole mass and 
hence the threshold is not tied to twice the pole mass onc e we talk about accuracies of several 
hundred MeV. One can make this cancellation manifested by subtracting the dangerously large 
terms from the potential and by adding them back to the pole mass. The result is a modified 
mass renormalization prescription, called the 'potential subtracted mass', which can be related 
to a more conventional definition like the MS mass by a reasonably well-behaved series, but 
which is at the same time not unphysically far away from the threshold like the MS mass (at the 
scale mt) itself. It is of course possible to use other modified mass renormalization prescriptions, 
provided they also lead to a manifest cancellation of long-distance sensitive terms and can be 
computed to sufficient accuracy. One possible alternative is discussed in Ref. e3. 

The discussion of the previous paragraph applies to a stable or unstable quark. Also, as 
mentioned, the top quark pole mass suffers from the same long-distance sensitivity as the bottom 
or charm quark pole mass despite the fact that the width of the top quark is significantly larger 
than Aqcd- Still, the width helps. The point is not that the top quark pole mass should be 
better behaved. The point is that contrary to bottom quarks, where one can find observables 
(such as the B meson mass) which are as long-distance sensitive as the pole mass, there is no 
observable involving top quarks that would be as sensitive to long distances as the top quark 
pole mass. In this precise sense, the top quark pole mass is an irrelevant quantity. The top 
quark propagator is simply never on-shell. The propagator is given by 1/(|( — rrit + iTt/2), but 
the loop or external momentum p can always be considered as real. The denominator of the 
propagator never gets smaller than about rritTt- We then find, for a quantity that would have a 
long-distance correction of (relative) order AqcB/fn for a stable quark. 



Aqcd _^ Aqcd _ Aqcd 
m m T 



(15) 



We expect, at least, a supDression by AQco/rj due to the finite width. (For a related discussion 
of this point, see Refs.EjcJ.) 
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Order 


n 


pole 




Ref. 


46 


NLO 


8 - 20 


4.827 ±0.007 


— 


Ref. 


47 




8 - 20 


4.604 ±0.014 


4.13 ±0.06 


Ref. 


48 


NLO 


10 - 20 


4.75 ±0.04 


— 


Ref. 


28 


NNLO 


10 - 20 


4.78 ± 0.04 




Ref. 


29 


NNLO 


4 - 10 


4.78 - 4.98 


4.16 - 4.33 


Ref. 


30 


NNLO 


14 - 18 




4.20 ±0.1 



Table 1: Bottom quark mass obtained from T sum rules in GeV. The MS mass is renormalized at the scale 
of its own value. The different values are not strictly comparable, because different values of as{mz) 
may be implied. Whenever available, we quote a result obtained from a fixed input rather than fitted 
as{mz). 'a^' means that no systematic resummation has been attempted. 



6 Bottom quark mass from Upsilon sum rules 

The derivatives of the bottom vector current two-point function at q'^ 
inclusive bottom cross section by 



127r2 



nl 



in 



ds 



are related to the 



(16) 



up to a (small) correction due to bb radiation from light quarks. The left hand side can be 
computed in perturbation theory; the right hand side from data. 

The parameters of the lowest T{nS) resonances are well- measured, but the continuum cross 
section near the threshold is not. Hence the experimental error of the right hand side decreases 
with increasing n, because higher moments weight lower s. But larger n makes the theoretical 
calculation more difficult, because the expansion parameter is Ogy/n- The enhancement by a 
factor y/n is a consequence of the Coulomb enhanced terms in the theoretical calculation of R. 
A power of f in i? corresponds to a power of n~^/^ in the moment. When 

\fTvnas{m\,l ^) <^ 1 (17) 

is no longer satisfied, a resummation by the methods discussed earlier in this talk is necessary. 
This is already the case for n larger than 4. There is an intermediate n where experimental and 
theoretical uncertainties are balanced. 

Recently, there have been several calculations which implemented this resummation at 
NNLO. Their results, together with results of previous (NLO or fixed-order a^) calculationsE§~@ 
are compiled in Table |l[ Note that as the accuracy of the calculation increases, the uncertainty 
in the result becomes larger. Instead of our own comments, we refer to Ref. cj for a discussion 
of this intriguing point. The result of Ref. is low, mainly because they do not include the 
sub-threshold poles of the Coulomb Green function into their calculation. But since to a large 
extent it is the average over Coulomb poles which is dual to the average of over the physical 
T resonances, their inclusion is necessary for moments which receive their largest contribution 
from the T resonances. Technically, this follows from the fact that the Coulomb poles contribute 
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and amount of order as{mh/^/n)^ to the moments, to be compared with the tree graph 

contribution. 

Note that because of the cancellation of long-distance contributions between the potential 
and the pole mass discussed in the previous section, it is advantageous to use a mass renormal- 
ization convention different from the pole mass beyond NLO. This strategy has been adopted 
in Refs.E3H. 

We mentioned earlier that there is a non-perturbative contribution to the heavy quark cross 
section at NNLO when mv'^ ~ Aqcd- There is a contribution to the moments from the scale 
nih/n. If we require for a perturbative treatment that all scales are larger than 0.5 GeV, this 
limits n < 10. The use of larger n in most of the calculations quoted in Table |l] has been 
justified by the observation that even at n = 20 the gluon condensate contribution to the 
moments is very small. Hence one may be in the fortunate situation that the actual non- 
perturbative contribution is smaller than the power counting argument suggests. However, the 
condensate expansion converges parametrically only when rrib /n ^ Aqcd • If the expansion does 
not converge, the size of its first term is not very conclusive. (For the Coulomb energy levels, 
the contribution from dimension-6 operators has been considered in Ref.E^.) 

7 Conclusion 

The past year has seen significant progress in our understanding of perturbative quark-antiquark 
systems close to threshold, both methodical and in terms of explicit higher-order calculations. 
We have come nearer to the answer to the question how accurately the bottom and the top 
quark mass can be determined by purely perturbative means. 
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